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We present realistic multiband calculations of scanning tunneling spectra in Bi2Sr2CaCu2 0g+« 
over a wide doping range. Our modeling incorporates effects of a competing pseudogap and pairing 
gap as well as effects of strong electronic correlations, which are included by introducing self-energy 
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I. INTRODUCTION 

A curious, topological feature of any two dimensional 
electronic band is a saddle-point Van Hove singularity 
(VHS), which can lead to a crossover from electron- like 
to hole like Fermi surfaces with doping, accompanied by 
a logarithmically diverging density-of-states (DOS) and 
a vanishing Fermi velocity^. While the divergence itself 
can easily be smoothened by disorder or interlayer cou- 
pling, the presence of a peak in the DOS is a robust fea- 
ture of the underlying spectrum. In hole-doped cuprates, 
most theoretical studies of doping dependence find en- 
hanced tendencies for superconductivity or competing 
phases - including the possibility of phase separation - 
when this VHS approaches the Fermi level Ep£& This 
scenario seems to work for La2- a; Sr a ;Cu04, where the 
VHS nearly coincides with Ep at optimal doping^, but in 
most other cuprates, with a larger second-neighbor hop- 
ping parameter t', optimal doping and the VHS doping 
seem to separate further as T c increases.— Notably, early 
angle-resolved photoemission spectroscopy (ARPES) ex- 
periments found evidence for a peak-dip-hump structure, 
where an electronic mode pushed a narrow coherent band 
with a prominent VHS ('peak') very close to the Fermi 
level of most cuprates, while pushing incoherent spectral 
weight away from Ep ('hump'), thereby leaving a 'dip' 
in the spectrum near the mode energy^ However, the co- 
herent VHS peak has proved elusive in these materials, 
and there is controversy as to whether it can be seen in 
tunneling spectroscopy 7 -^. 

Here we explore this issue by comparing recent scan- 
ning tunneling microscopy and spectroscopy (STM/STS) 
data with a realistic multiband modeling of the spectra 
in Bi2Sr 2 CaCu20 8+ «5 (Bi-2212). The model is able to re- 



produce the STM features in detail, shedding light on the 
role of the VHS and the ordering phenomena involved in 
'two gap' physics. 

Our doping-dependent tunneling spectrum is calcu- 
lated using a multiband Nambu-Gorkov Green's func- 
tion tensor formalism. The main technical innovation 
of this work is to incorporate a variety of intermediate 
coupling effects into the tunneling formalism via a tensor 
self-energy extracted from quasiparticle-GW (QP-GW) 
calculations^. We model the two competing order gaps 
as due to BCS-type SC couplin g 10 ' 11 and antiferromag- 
netic order. To these we add intermediate coupling self- 
energies, calculated in a GW model. Our QP-GW self- 
energies have been shown to capture a wide range of phe- 
nomena in cuprates over an extended range of doping and 
excitation energies, including the relationship between 
magnetic order and superconductivity, and the so-called 
waterfall physics*^ Notably, for a realistic consideration 
of the features in STM/STS, it is essential to utilize a 
multiband framework. This makes it possible to incor- 
porate filtering and tunneling path effects due to the BiO 
and SrO layers separating the vacuum and the cuprate 
layers, and the important role of the apical oxygen and 
Cu d z 2 orbitals^^i 

While the one band QP-GW model employed to obtain 
the self-energies used in this study is nearly parameter- 
free^, the present calculations unavoidably introduce ad- 
ditional parameters due to two separate factors. Firstly, 
the pairing gaps are treated as parameters, since the 
Eliashberg theory cannot accurately predict the gap 
size, particularly in the pseudogap phase. A second set 
of parameters arises because STM experiments on the 
cuprates find that the gap size varies in different local 
patches across a given crystal. In our modeling, we as- 
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sume that different patches represent regions of different 
local doping^—. A strong inference from our study is 
support for this local doping model since a number of 
observed spectral features can be correlated with theory 
over a wide range of energies and dopings. 

This article is organized as follows. Section II de- 
scribes our model of electronic spectrum and gives a brief 
overview of the formalism of tunneling calculations. Sec- 
tion III presents theoretical results for the normal and the 
superconducting state as a function of doping and com- 
pares and contrasts these results with the correspond- 
ing experimental results, including issues related to the 
two-gap physics and signatures of the VHS. Section IV 
discusses broader implications of our study. This is fol- 
lowed by a summary of our conclusions in Section V. 
Appendix A gives details of the interaction terms used 
in the Hamiltonian, while Appendix B discusses details 
of the self-energy corrections used for incorporating elec- 
tron correlation effects in our tunneling computations. 
Appendix C summarizes a related model which displays 
peak-dip-hump physics. 



II. METHODOLOGY 

A. Green's function formulation of tunneling 
current 

Our computations are based on a multiband tight- 
binding model using 4 x 116 orbitals per unit cell, i.e., 
a fourfold basis taking into account electrons and holes 
for both spin directions! 10 ' 11 For treating the magnetic 
order and superconducting pairing, the tensor (Nambu- 
Gorkov) Green's function Q is employed with the corre- 
sponding Dyson's equation: 19 
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where G e and Gh, denote the Green's functions for the 
electrons and holes, respectively 2 ^, and the matrix ele- 
ments of operator V represent the interaction terms of 
Hamiltonian of Eq. (7), below. 

In the tunneling calculations the filtering effect of the 
surface layers is taken into account by using the Todorov- 
Pendry approach 2 ^ in which the differential conductance 
a between orbitals of the tip (t, t') and the sample (s, s') 
is given b y 10 ' 11 
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is expressed in terms of the retarded and advanced elec- 
tron Green function or propagators. Eq. [5] differs from 
the more commonly used Tersoff-Hamann formulation 22 
in that it takes into account the details of the symme- 
try of the tip orbitals and how these orbitals couple with 
the orbitals of the cuprate layer through the filtering BiO 
and SrO layers. 

Equation [3] can be rewritten as^i 



27re 2 x - 

tt'cc' 



where 



ptt, [E F )M VcPcd (E F + eV)Ml v (4) 



M tc = V ts G°+V fc , (5) 



gives the filtering amplitude between the cuprate layer 
and the tip, and constitutes a multiband generalization 
of the filtering function of Ref. HI- Note that the matrix 
element of the density of states operator p CC ' within the 
cuprate plane can be recovered in terms of the spectral 
function: 

Pec = -- E( G -£"«G- C , + F+V' a F- c ,) (6) 



Notably, the matrix M itself is rather structureless, but 
it is very selective in the way it couples the insulating 
layers and the cuprate layer due to the term V/ c ~> In 
fact, for symmetry reasons the tunneling channel passes 
through the d z i orbital and not the d x 2_ y 2 orbital of the 
Cu atoms. 



B. Interaction Effects 

In order to meaningfully describe the doping depen- 
dence of the tunneling spectra of the cuprates, various 
key interactions must be included in the Hamiltonian. 
To this end, we introduce interactions via a Hubbard U 
on the Cu sites, which are approximated by a mean-field 
antiferromagnetic order plus GW-type self-energy asso- 
ciated with spin and charge fluctuations. We model the 
two competing order gaps as being due to a BCS-type SC 
coupling, and an antiferromagnetic order. Superconduc- 
tivity is included by adding a d-wave pairing interaction 
term A and the AFM order— is included by introducing 
a Hubbard term Um, leading to the Hamiltonian: 



(2) 



with the real-space creation (annihilation) operators c' a(T 
(c Q(J ). Here a is a composite index for the type of orbital 
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FIG. 1: (Color online) Comparison of the theoretical (solid) 
and experimental (dashed, Ref. l28l ) normal state STS spectra 
over the hole doping range x — 0.21 — 0.17. 



and its site, and a is the spin index. e a is the on-site 
orbital energy, while V a p is the hopping integral between 
orbitals a and f325r— The band structure related param- 
eters e a and V a p are the same as those used in Ref. [Tl|. 
Details of the interaction parameters U, m, and A are 
given in Appendix A. 

In order to model effects of the low and high energy 
bosonic couplings and background scattering in the elec- 
tronic system, we add a self-energy correction in the 
propagators, which consists of three distinct contribu- 
tions: 

S Q ^( £ ) = E^(e) + E^(e) + ^J{e). (8) 

Here, superscripts L and H refer to corrections due to 
the low- and high-energy bosonic couplings, respectively, 
while the superscript imp refers to the correction term 
arising from background impurity scattering. For S L 
we use a phonon contribution based on a Debye spec- 
trum and Yi H is obtained within the framework of a self- 
consistent GW scheme. Details of the self-energy cor- 
rections of Eq. 8, including an overview of the GW- 
schcme used, are given in Appendix B. We note that 
in the theoretical STS spectra discussed in Section III 
below, the electronic contribution T, H dominates in pro- 
ducing key features in the spectra, while £ im P yields a 
realistic broadening of the spectra and S L contributes to 
the sharpness of the coherence peaks. 



III. RESULTS 
A. Normal State Spectra 

Fig. 1 compares the normal state spectra taken above 
the superconducting transition temperature T c where su- 
perconductivity is absent although a pseudogap could be 
present. Experimental data over a wide energy range are 
available only near and above optimal doping where the 
pseudogap is small or zero.— In Ref. HH, the supercon- 
ducting gap map of an inhomogenous sample was first 
measured at low temperature to determine the local gap 
A in different domains on the surface. The temperature 
was then raised above the highest superconducting tran- 
sition temperature, and the corresponding normal state 
spectra were measured. Along these line o 14 ' 15 , we use 
the measured low-T gap values varying from 20 meV to 
36 meV to adduce the local doping, which yields dop- 
ing values varying over x = 0.21 — 0.17 as shown in Fig. 
1. The theoretical as well as the experimental spectra 
in Fig. 1 are seen to be quite featureless, except for the 
pronounced asymmetry between high positive and nega- 
tive bias voltages, which we have analyzed previously^ 
The most prominent feature in the spectra is a large 
hump-like feature located around -200 meV binding en- 
ergy for x= 0.21, which shifts away from the Fermi level 
to higher binding energies and becomes smoother as the 
doping is decreased. Our calculations are seen to repro- 
duce this hump-like feature and its characteristic shift 
and smoothening with doping in substantial accord with 
experimental results. It should be noted that over the 
doping range considered in Fig. 1, pseudogap effects in 
our computations are quite small. We have verified this 
by carrying out a series of computations where the pseu- 
dogap was artificially turned off (not shown in Fig. 1 for 
brevity), suggesting that the pseudogaps have closed at 
the temperature of 93K at which the experimental data 
are taken. 

Figure [^a) gives a plot of the energies of various fea- 
tures in the theoretical and experimental STS spectra 
of Fig. 1 to allow a more quantitative discussion of the 
doping dependencies of these features and their correla- 
tion with the phase diagram of Bi2212. The model in- 
cludes bilayer splitting, and the solid blue [short-dashed 
red] curve represents the VHS peak of the antibonding 
[bonding] band. We see that the experimental normal 
state hump (blue filled circles)^ is in good agreement 
with the doping dependence of the antibonding VHS. The 
calculated VHS position does not change significantly in 
the superconducting state, so we also plot the experimen- 
tal hump feature from the superconducting state (violet 
filled diamonds) 17 , from Fig. 3 below. When an AF pseu- 
dogap turns on, a second feature appears [green long- 
dashed and orange dotted lines for the antibonding and 
bonding bands respectively]. Frames (b) and (c) explain 
the origin of this feature for the antibonding band, la- 
beled B. These figures are based on a three-band model, 
containing a Cu d x 2_ y 2 and two O p orbitals, in which 
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FIG. 2: (Color online) (a) Derived AF state phase dia- 
gram. The filled blue circles [Ref. [2S| and violet diamonds 
[Ref. Q3] are the hump features for Bi2212 derived from re- 
cent STM measurements. These are compared to a calculated 
doping dependence of the bonding [B] and antibonding [A] 
VHS (red short-dashed and blue solid lines, respectively) in 
Bi2212, based on the present analysis, (b) Calculated dis- 
persion at x = 0.15, with feature A, the VHS of the lower 
magnetic band, corresponding to the blue solid line [or red 
short-dashed line, for the bonding band] in (a), and feature 
B, the bottom of the upper magnetic band, corresponding to 
the green long-dashed line [orange dotted line] in (a). The 
thickness of the lines represents the spectral weight due to 
the AF structure factor, (c) Corresponding DOS, showing 
features derived from A, B, C, and D in (b). (d) Blow-up of 
the near-i?F dispersion. 



bilayer splitting is neglected, to clarify the origin of the 
features. From Fig. 3(b), it can be seen that the AF gap 
has a strong fc-dependence, splitting the band into upper 
and lower magnetic bands (U/LMBs). Features A and D 
are the VHSs of the LMB and UMB respectively while 
feature B is the bottom of the UMB and feature C the 
top of the LMB. Note the characteristic form of the DOS 
associated with each feature in Fig. 2(c). 

As the AF pseudogap shrinks to zero, features A and 
B [and C and D] merge. That is why there is no near-_Ep 
feature in the three spectra in Fig. 1 with higher dopings. 
For lower doping a feature corresponding to B appears 
in the low-T theoretical curves, but seems to be absent 
in the T=93K experimental data of Fig. 1, suggesting 
a T-dependent pseudogap closing. Note that at lower 
doping feature B moves closer to the Fermi level, and 
should persist to higher temperatures, as the pseudogap 
is larger. 



B. Superconducting State spectra 

Figure |3ta) compares our theoretical results^ with 
the corresponding experimental spectr a ' 30 at five rep- 
resentative values of hole doping, with doping decreasing 
from the optimally doped SC phase at x — 0.19 at the 
top to the underdoped pseudogap phase for x — 0.08 at 
the bottom of the figure. With reference to Fig. 1, we 
see that the hump-like feature at negative energies and its 
characteristic doping dependency in the normal state is 
present also in the superconducting state spectra of Fig. [3] 
(plotted in Fig. [21a)). The main difference between the 
normal and superconducting state spectra is the opening 
up of a complex gap structure at low energies, which is 
highlighted in the expanded view of Fig. \T[b) We note 
here with reference to the lowest set of curves in Fig. 3(a) 
as an example that the STM spectra resemble the antin- 
odal ARPES spectra.— Moreover, the angle-dependence 
of the gap function adduced from the STM quasiparti- 
cle interference maps is in good agreement with the gaps 
obtained via ARPES measurements. 33 These results sup- 
port our use of the local gap to determine the local doping 
of the patches in STM data. 

Figure 2] gives a plot of the two low-energy gaps 
from the theoretical STS spectra of Fig. 3, shown as 
filled and empty blue diamonds. These are seen to 
be in good agreement with the 'two gap' features ex- 
tracted from ARPES^^ and with Raman results^&^I for 
HgBa 2 Cu0 4+5 (Hgl201). The remaining curves 3 ^ will 
be discussed below. Interestingly, our computations re- 
produce the remarkable anticorrelation displayed by the 
large hump feature at negative energies in the experi- 
mental STS spectra in Figs. 2 and 3 in the way it varies 
with doping in relation to the doping dependence of the 
size of the gap2£. Specifically, with increasing doping, 
as the hump moves closer to the Fermi energy and the 
density of states at the Fermi energy increases, the size 
of the gap in Figs. 3 and 4 does not become larger as we 
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FIG. 3: (Color online) (a) Theoretical (solid) and experi- 
mental (dashed, Ref. 117]) superconducting state spectra over 
the hole doping range from x = 0.19 at the top to x = 0.08 
at the bottom of the figure.— 9 For comparison, the experi- 
mental ARPES spectrum^ at x = 0.08 (black dashed curve) 
is shown. The dot-dashed line is drawn through the hump 
feature in the computed spectra as a to guide the eye. (b) 
Blowup of theoretical spectra near the Fermi level is com- 
pared with the corresponding data of Ref. |3Cj. The black and 
red bars indicate the gap widths shown as empty and filled 
diamonds in Fig. 4. 



would expect from the BCS theory, but instead the gap 
becomes smaller. The size of the coherence peaks, on the 
other hand, is seen to grow with increasing doping in our 
theory in accord with experiments. 

Fig. 5 gives insight into the origin of the two gap 
physics in Bi2212. For this purpose panel (a) gives a 
map of the spectral weight obtained from the imagi- 
nary part of the one-particle Green's function along the 
r — X(ir, 0) — M(tt, 7t) — T line in the k-space. This map 
may be thought of as an effective theoretical E vs k spec- 
trum resulting after the self-energy correction £ has been 
applied to the propagators. In the absence of self-energy 
corrections, this map will consist of 5-function peaks aris- 
ing from the real poles of the Green's function corre- 
sponding to the familiar energy bands of Bi2212 modi- 
fied by the the AFM and SC orders .-^ii, At the doping of 
x = 0.08 considered in Fig. 5, the antibonding VHS lies 
around -300 meV, well below the Fermi energy. After self- 
energy corrections are introduced, the VHS is broadened 
substantially and yields the large hump-like structure in 
the local density of states around -300 meV as seen from 
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FIG. 4: (Color online) Two gap physics in the superconduct- 
ing state. The two theoretical gaps from Fig. 3(b) are plotted 
as filled and empty blue diamonds and compared to exper- 
imental ARPES results from Ref. [34] (green triangles) and 
Ref. dH (brown inverted triangles) and corresponding Raman 
results for HgBa2Cu04+5 (Hgl201) (black long- and short- 
dashed lines) from Ref. [37] . Shown also is the input SC gap 
parameter, violet squares, the expected position of a VHS, 
based on an experimental ARPES dispersion- orange dotted 
line, and a calculated SC red circles, to illustrate the 

doping dependence expected for a spin-pairing mechanism. 



the LDOS plot of Fig. 5(b) [leftmost arrow]. Due to bi- 
layer splitting 42-46 , the bonding VHS lies near -800 meV 
at this doping, and is not seen in Fig. 5. 

Next, we address the features near Ep in Figs. 3(b) 
and 4 and how these features relate to the two- 
gap physics, which has been discussed in connection 
with Raman scattering^^ 7 , angle-resolved photoemis- 
sion (ARPES)^ 4 .^, and recent STM studies^L^. Two 
different components of the gap have been observed with 
different doping dependences^: a nodal pairing gap A„ 
with a parabolic doping dependence reminiscent of the 
superconducting T c , and an antinodal gap A* which 
increases roughly linearly with underdoping, similar to 
the pseudogap onset temperature T*. The experimen- 
tal near-gap features^ in Figs, [^b) and 4 clearly reflect 
the presence of two gaps, which are well captured by the 
present calculations. While the overdoped sample dis- 
plays a single gap feature, samples at lower doping show 
four features comprising two gaps, an outer gap which 
grows with underdoping and an inner gap which shrinks 
(see Fig. 4), in good agreement with the two-gap sce- 
nario. Coherence peaks are prominent in the overdoped 
sample, but gradually disappear as doping is lowered, in 
good accord with the computed STM spectra. [In the un- 
doped case we find a large magnetic gap, with no evidence 
of superconductivity (not shown).] From Fig. 3(b) it can 
be seen that our modeling reproduces the electron-hole 
asymmetry of the experimental coherence peaks, which 
is unexpected in the BCS formalism^. 

Fig. 6 gives the decomposition of the LDOS of Fig. 
5(b) for x=0.08 into the spin-resolved contributions from 
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FIG. 5: (Color online) (a) Calculated E — k spectral weight 
diagram for x = 0.08 along high symmetry line F — (n, n) — 
(-7T, 0) — r in momentum space. White/brown color denotes 
high and blue color low intensity, (b) The corresponding total 
LDOS from Cu sites, with arrows relating three prominent 
LDOS peaks to spectral features in (a). Note that several 
peaks seen in LDOS in the low energy region do not show up 
in (a) because the associated spectral weight lies along the 
zone diagonal, which is not shown in (a) [see Fig. 6(b)]. 



the upper and lower magnetic bands (UMB and LMB). 
In our computations there is a small magnetic gap sep- 
arating the UMB and LMB, leading to electron pockets 
near (tt,0) and hole pockets near (■§,§)• The two gaps 
in the total LDOS curve in Fig. 6 are identified by the 
arrows labeled 1 and 2, and are plotted in Fig. @] as blue 
diamonds. In order to clarify the role of pairing in these 
two features, Fig. [B] also shows the anomalous matrix 
elements F a p of the Nambu-Gorkov Green's function 51 , 
which can be used to monitor the pairing effect. The con- 
tribution of the anomalous elements to both the in-gap 
and gap-edge peaks indicates that pairing is involved in 
both features. Notably, the gap features in Fig. 3(b) more 
closely resemble the UMB features in Fig. 6, indicating 
that the tunneling matrix element couples differently to 
the LMB and UMB states. 

Further insight into pairing can be obtained by con- 
sidering the evolution of dispersions with binding en- 
ergy. Fig. 7 shows spectral maps in momentum space 
at the energies of the coherence peaks for the lower edge 
of the pseudogap in frame (a) and for the lower in-gap 
peak in frame (b). The corresponding anomalous spectral 
weights are shown in frames (c) and (d). The peak at the 
pseudogap edge is associated with the top of the (tt, 0) 
pocket as can be seen from the strong pairing intensity in 
Fig. 7(c). The in-gap peak is associated with the pairing 
gap on the (ir/2, n/2) pocket, Fig. 6(d), where AF order 
is absent. Note that at this energy the (it, 0) pocket is ab- 
sent, Fig. 7(b). The (7r, 0) peak is found to predominantly 
represent AF order with a Bilbro-McMillan-like dressing 
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FIG. 6: (Color online) Calculated spin resolved LDOS of the 
Cn-d x2 -y2 orbital for x = 0.08: Total LDOS (blue solid line) 
and its decomposition into contributions from the lower mag- 
netic band (LMB, red dashed line) and the upper magnetic 
band (UMB, green dashed line). Also shown are the corre- 
sponding imaginary parts (cyan and magenta solid lines) of 
the anomalous matrix elements F a p of the Nambu-Gorkov 
Green's function for the d X 2- y 2 orbitals on two neighbouring 
Cu atoms. This matrix element indicates strong electron-hole 
hybridization or pairing effect at the pseudogap edge labeled 
by '1' and the in-gap peak labeled by '2' in the total LDOS 
(solid blue line). 



by superconductivity 5 ^, while the (tt/2,tt/2) peak is a 
pure SC gap. Note in Fig. 7(b) that the spectral weight 
of the (7r/2,7r/2) pocket resembles an arc with peak in- 
tensity at the AF zone boundary as seen in STM stud- 
ies. With decreasing doping the gap between the LMB 
and UMB increases, leading to a monotonically increas- 
ing pseudogap energy. In contrast, the Fermi surface on 
the LMB near (ir/2, tt/2) shrinks monotonically since its 
area is proportional to x. For a <i-wave gap which van- 
ishes at (7r/2, 7r/2), its maximum value must ultimately 
decrease with underdoping, thereby explaining the oppo- 
site doping dependence of the two gaps. 



C. Comparison with Bi2201 

To gain more insight into the relation of Eyns and 
the superconducting gap A, in Fig. 8 we compare the 
present results with recent experiments on Bi2Sr2Cu06+<5 
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FIG. 7: (Color online) (a) Spectral map at constant energy 
for x = 0.08 for the regular matrix elements of the Nambu- 
Gorkov Green's function. The energy corresponds to the lower 
edge of the pseudogap, i.e. the energy of the peak marked 1 
in Fig. 5. (b) Same as (a), except the energy corresponds to 
the energy of peak 2 in Fig. 5. (c) Same as (a), except it 
gives the spectral map for anomalous matrix elements of the 
Green's function at the energy of peak 1. (d) Same as (c) for 
the peak 2. Arcs in (c) and (d) have zero weight along the 
nodal line, where the gap vanishes. 



The decrease of A with increasing Evhs, Fig. 8(b), is 
quite similar in both materials, although with somewhat 
faster change in Bi2201 (note change of scale). The re- 
markable resemblance to the dotted curve is unexpected, 
and deserves some comment. The dotted curve repre- 
sents the relation 

A 2 + A 2 = A\E VHS \, (9) 

with Ao ~ 18 meV, and A ~ Aq/2. The general form 
Y) . Af — constant is a common finding for competing 
orders, but in Eq. [5] the meaning of Ao is unknown, and 
the right-hand side is not constant, but vanishes when 
the VHS is at the Fermi level. This is suggestive of yet 
another quantum critical point (QCP) in the cuprates, 
but this one associated with the VHS. 3 ' 54 Unlike most 
VHS-related features, in this case the gap scales to zero 
at the VHS, suggestive of nesting of the antinodal parts 
of the Fermi surfaced 

Remarkably, QCPs have been proposed at three dop- 
ings in the cuprates: in the underdoped regime, where 
superconductivity is lost in a superconductor-insulator 
transition, near optimal doping, and now in the over- 
doped regime, where superconductivity vanishes near a 
VHS. Equation [9j or its generalization to multiple order 
parameters, offers a rationale for such a series of transi- 
tions, as the shrinking nesting inaugurates an ever more 
desparate competition among phases. 

IV. DISCUSSION 
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FIG. 8: (Color online) (a) Tunneling dl/dV spectra for 
Bi2201 corresponding to Evhs = -10 meV (upper, red line) 
and -56 meV (lower, violet line, offset for clarity), after Ref.[8|. 
(b) Plots of A vs Evhs for Bi2212 (red diamonds and vio- 
let triangles from Ref. and blue circles from Ref. l28h and 
Bi2201 (inverted green triangles from Ref. [H). Dotted line 
indicates a square-root dependence of A on Evhs + 36 meV. 



(Bi2201)4 We see that the shape of the VHS at highest 
doping is consistent with the calculated shape in Bi2212 
at lower doping compare the top curves in Fig. 3(a) and 
Fig. 8(a) but that it gets considerably broadened for 
slightly lower doping, lower curve in Fig. 8(a). This may 
explain why our calculations in Fig. 3(a) underestimate 
the VHS broadening, as the QP-GW technique tends to 
produce too little broadening at higher energies. 



The present calculations strongly support the idea that 
a competing order pseudogap provides a good description 
of two gap physics in hole doped cuprates The model 
offers new insights into the role of the VHS. It has been 
widely debated whether tunneling measurements can see 
the VHS°££. The results are quite revealing and show 
that over the entire doping range studied, the feature 
conventionally identified as the VHS (orange dotted line 
in Fig. 4) is actually a gap feature^ 5 " 5 " which can be traced 
to the bottom of the upper magnetic band. This explains 
why no VHS is found in normal state tunneling near op- 
timal doping when the superconductivity is turned off. 
The calculation also provides insight into the transition 
to a large-gap insulator at half filling. The AF pseudogap 
is essentially given by the separation between the VHS 
and the bottom of the upper magnetic band, features A 
and B in Figs. 2(b), 2(c). From Fig. 2(a), this can be 
seen to be approximately 0.3 eV at x = 0.08, and the 
route to a 2 eV gap at half filling is clear. 

Instead, we find the true VHS feature at higher ener- 
gies. The doping evolution is consistent with that found 
in overdoped Bi2201 cuprate^ where the VHS is clearly 
seen close to the Fermi level. The anticorrelation of the 
VHS peak with superconductivity could arise because the 
VHS can drive a competing ferromagnetic instabilit y 39 ^ 6 

It should be noted that the present model is not con- 
sistent with the conventional picture of peak-dip- hump 6 . 
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By a relatively small change of parameters we have found 
a second solution, much closer to this conventional form. 
While this model reproduced many of the features of the 
STM data, it was ultimately less satisfactory than the 
model described in the text. In Appendix C, we briefly 
review this peak-dip-hump model. Among other prob- 
lems noted in Appendix C, the peak-dip- hump model 
is unable to explain the transition to a large-gap insu- 
lator at half filling, producing a magnetic gap of only 
~100 meV even at the lowest doping. 

While cuprates have traditionally been treated in 
strong coupling formalism, recent investigations have 
shown that intermediate coupling models, such as the 
present QP-GW model, can capture many salient fea- 
tures of cuprate physics, including the doping depen- 
dence of dispersion and optical properties^ 7 - - — . In par- 
ticular, our QP-GW self-energy successfully describes 
the dressing of low-energy quasiparticles by spin and 
charge fluctuations,— the high-energy kink (HEK) seen 
in ARPES, 60 the residual Mott bands in the opti- 
cal spectra^ gossamer features^ anomalous spectral 
weight transfer with doping^, and the magnetic reso- 
nance in neutron scattering^ In fact, our intermediate 
coupling model yields a number of characteristic signa- 
tures of strong coupling physics of the Hubbard model, 
including suppression of double occupancy^, the t—J- 
model dispersions^, spin wave dispersion™,, the 1/U 
scaling of the magnetic order^l 7 - and the phenomenon 
of anomalous spectral weight transfer (ASWT) 62 i 68 ' 69 . 
Interestingly, a recent variational calculation finds a 
smooth evolution from a SDW to a Mott gap with in- 
creasing C7, with no intervening phase transition or spin 
liquid phase in the cuprate parameter ranged 



V. SUMMARY AND CONCLUSIONS 

We have carried out a realistic, multiband modeling of 
STS spectra in Bi2212 over a wide doping range where ef- 
fects of competing pseudogap and superconducting gaps 
are incorporated, and strong electron correlations are 
treated by introducing self-energy corrections in the one- 
particle propagators. The theoretical spectra capture 
many of the salient features of the corresponding exper- 
imental spectra and their doping dependencies, includ- 
ing the two-gap behavior at low energies, electron-hole 
asymmetry of the coherence peaks, and the presence of a 
prominent hump feature at high energies. Our analysis 
yields insight into the complex manner in which the VHS 
manifests itself in the STS spectra. 
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Appendix A: Details of interaction terms in the 
Hamiltonian 

Here we provide details concerning the superconduct- 
ing and antiferromagnetic terms in Eq. [7] The mean field 
superconducting coupling between electrons and holes is 
of the form 

A Q( 3 = yU a p a b{Cg\.Cb\) ■ (Al) 
ah 

We assume a (i-wave form for the gap as is appropriate 
in the cuprates.— The interaction U a p a b is not known, 
and hence we follow the common practice of introducing 
a gap parameter, which gives the correct gap width and 
symmetry 7 ^. The gap values of | A| = 25 — 45meV, plot- 
ted as violet squares in Fig, 4, are chosen to model the 
experimental spectr It is interesting to note that 

the gap calculated from a spin fluctuation mechanism 39 
(red circles in Fig. 4) has a very similar doping depen- 
dece, andeed a very similar magnitude to the smaller gap 
feature. 

Based on the QP-GW results, we incorporate magnetic 
interactions into our model via an onsite interaction term 

H = U aa n a ^n a ^ = U aa c^ a ^c a ^c^ a ^c a ^. (A2) 

At mean field level, this reduces to a magnetic interaction 

where m a = n a ^ — n a ^. In the present calculation, 
we assume a doping independent U = 7.5 eV and the 
< m a > are found self-consistently from a mean field 
calculation ! 29 ' 44 ' 72 



Appendix B: Details of Self-energy corrections 

Details of the three contributions to the self-energy 
of Eq. [8] are as follows. In the phonon contribution E L , 
the phonons are approximated by a Debye spectru m 73 ' 74 . 
The self-energy can then be written as a convolution of 
the Debye spectrum with a constant density of states for 
the electrons as^ 

E*(e) = ((2z + in) + (z 2 - l) In ) , 

(Bl) 

where z = (e + irj) / (hCld) , A = ^-T 2 p, and rj is a con- 
vergence parameter. We assume the self-energy for sim- 
plicity to be of a diagonal form, and apply this to Cu- 
d x 2_ y 2 orbitals using h£ld = 80meV and A — 60meV. 
When combined with these values yield a reason- 
able description of the peak-dip-hump structure and the 
smoothness of the spectrum. 
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For the electron-electron self-energy, X , we utilize a 
simple Fermi-liquid type self-energyiS, except for the an- 
tibonding C11O2 band, nearest the Fermi level, for which 
we introduce a QP-GW self-energy^ 

E ff (fc,c) = E/ ^(A-ftw + w'Wg,^) 

where W — (U 2 /2)Im(2>Xs + Xc) and U is appropriate to 
the one- band model. The RPA spin/charge susceptibil- 
ities are Xc/s(<?,w) = xo(g,w)[l ± f7xo(9, and the 
bare two-particle correlation function xo(<7, w) is the con- 
volution of the green's function. Self-consistency in X is 
achieved by calculating an average renormalization fac- 
tor Z = {l-dY./du)- 1 with the QP-GW method^^^i. 
The vertex correction is T = 1/Z. 

For electron doped cuprates, the competing order is 
(7r, 7r) antiferromagnetism, while the origin of the pseu- 
dogap is unclear in the hole-doped cuprates. We have 
shown- that the tunneling spectra are insensitive to the 
exact nature of the pseudogap, as long as it represents a 
density wave order [charge, spin, or flux phase] at (7r,7r), 
so here we model the pseudogap as a (tt, it) AFM. In- 
cluding a d-wave SC gap below T c , the green's function 
becomes a 4 x 4— tensor. The QP-GW self-energies are 
calculated in a (|fc), \k + Q)) basis, then transformed via a 
unitary transformation to a real-space basis involving the 
spin up and spin down sublattices. In this transformed 
self-energy matrix, the off diagonal terms are found to be 
small and of varying sign, and are neglected. The corre- 
sponding matrix elements for the spin up and spin down 
holes are obtained by using the relations derived in Ref. 

u 

Finally, an impurity self energy S jmp (e) = —i5.0meV 
is used in the SC phase for all the orbitals, which 
produces a reasonable smearing in the low-energy re- 
gion. For the normal state, a larger value of E"™ p (e) = 
— i20.0meV is used to account for enhanced normal state 
broadening near the Fermi energy. 

Appendix C: Peak-Dip-Hump Model 

The calculations described in the main text involve 
a moderate band renormalization by correlation effects, 
corresponding to a renormalization factor Z = 0.5, sim- 
ilar to that found in LSCO and NCCO. However, in 
Bi2212 ARPES experiments suggested a much narrower 
coherent band, with Z ^0.28^- In this Appendix we 
briefly describe the results of a model calculation with 
this smaller Z value, which also reproduces many of the 
features of the data of Figs. 1-4. 



This small-Z model reproduces a form of peak-dip- 
hump scenario. The strong renormalization means that 
the spectral weight below the Fermi level is split into a 
broad, weakly renormalized incoherent band (hump) and 
a narrow, strongly renormalized coherent band (peak). 
In this model, the broad feature at high energies, with 
the dispersion in Fig. 2, is the VHS of the incoherent 
band, while the coherent band has a VHS at much lower 
energies, which varies with doping exactly as the VHS 
derived from the ARPES dispersio n 38 ' 75 , orange-dotted 
line in Fig. 4. When an AF pseudogap and d-wave su- 
perconductivity are included, this model also reproduces 
the two-gap physics of Fig. 4. The dip between these two 
VHS features closely follows the measured position of the 
magnetic resonance peaki2&. 

There are several experimental features that this small- 
Z model does not successfully reproduce. For instance, 
the coherent VHS peak is always found to be more in- 
tense than the incoherent VHS, contrary to experiment, 
Fig. 3. This is not unexpected, as the QP-GW model is 
known to underestimate the incoherent spectral weight. 
However, there are additional problems, involving only 
the low energy, coherent spectral weight. Most impor- 
tantly, there is a true, coherent VHS near the Fermi level, 
which should clearly appear as a low energy peak when 
the AF and superconducting orders are turned off, in 
sharp contradiction to the normal state results plotted 
in Fig. 1. A more subtle feature involves the electron- 
hole asymmetry of the larger of the two-gap features as 
a function of doping. In this doping range, the coherent 
VHS always lies below the Fermi level, so the negative 
energy peak should always be more intense, and indeed 
the peak asymmetry should increase with underdoping, 
as the VHS moves further from the Fermi level. How- 
ever, at very low dopings, the experimental asymmetry 
reverses sign, Fig. 3(b). Indeed, this can be understood 
from the model calculations of Fig. 2(a). Feature B, the 
lower edge of the UMB, actually crosses the Fermi level, 
lying above it at low dopings. With this crossing, the 
theoretical asymetry reverses, in good agreement with 
experiment, Fig. 3(b). Finally, while the large-Z model 
clearly extrapolates to a large-gap insulator at half fill- 
ing, the small-gap model does not. The magnetic gap is 
directly reflected in the two gap physics near the Fermi 
level. Hence, even at the lowest doping it remains of or- 
der ^100 meV, and the manner in which a ^2 eV gap 
arises at half filling remains a puzzle. It is for these rea- 
sons that we prefer the larger- Z model. The apparently 
accidental pinning of the bottom of the UMB near the 
Fermi level may actually be a signature of a cooperative 
interaction between the AF and superconducting order, 
thereby stabilizing both phases. 
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